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We re-analyze Tavatron data on charmonium hadroproduction in the framework of the color-octet model 
implemented in the event generator PYTHIA taking into account initial-state radiation of gluons and Altarelli- 
Parisi evolution of final-state gluons fragmenting into cc pairs. We obtain new values for the color-octet matrix 
elements relevant to this production process. We discuss the sensitivity of our results to the transverse momentum 
lower cut-off employed in the generation to avoid the problematic pt—*0 region, arguing about the reliability of 
our previous extraction of the NRQCD matrix elements for the 3 S^ 8 ' and 1 5g 8 ' + 3 Pj 8 ^ contributions. Finally we 
extrapolate to LHC energies to get predictions on the J/ip direct production rate. 



1. INTRODUCTION 

In a series of previous papers we ana- 

lyzed charmonium hadroproduction in the light of 
the color-octet model ||. We employed a Monte 
Carlo generator (PYTHIA 5.7) § to assess the 
importance of some higher-order QCD effects. As 
a consequence we concluded that NRQCD ma- 
trix elements M determined from Tevatron exper- 
imental data B were considerably overestimated 
in other analysis H] , since they absorb some per- 
turbative effects, e.g. initial-state emission of glu- 
ons. The latter gives rise to an effective transverse 
momentum, enhancing the high-j> t tail of the dif- 
ferential distribution of charmonium. Once taken 
into account, the extracted long-distance param- 
eters have to be lowered significantly. This is es- 
pecially apparent for the linear combination of 
matrix elements (MEs) 1 Sq 8 ' 1 + 3 Pq 8 ^ where, in 
fact, a large discrepancy w.r.t. HERA results on 
J/ip photoproduction has been pointed out fT(| . 

However, in our former work we did not con- 
sider Altarelli-Parisi (AP) evolution of the frag- 
menting gluon into charmonium for the 3 S^ con- 
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tribution at high pt- In this paper we implement 
such an effect benefitting from the evolution per- 
formed by PYTHIA itself as will be explained be- 
low. We shall show that although our previous 
numerical values for the 3 s[ 8 ^ ME were somewhat 
low, no significant change occurs for the corre- 
sponding + 3 P (8) ME. We also discuss the 
sensitivity to the lower cut-off used to avoid the 
singular pt^O region. We get started with the 
latter point in the following section. 

2. SENSITIVITY TO THE p[ nin CUT- 
OFF 

Production cross sections which are singular at 
Pt = are automatically regularized in running 
PYTHIA by setting a p t lower cut-off at 1 GeV 
|(| . Consistently with our aim of using by-default 
options whenever possible in the generation, we 
kept this value in all our previous analysis Jl],|[ . 

In fact, the differential cross section can have 
a singular behaviour at low p t even to next-to- 
leading order due to imperfect cancellation of 
real and virtual contributions. However, from a 
NLO calculation of charmonium production by 
Mangano ]Tl[ |, one may conclude that our previ- 
ously used lower cut-off (p™ n = 1 GeV) is con- 



sistent with the threshold obtained by Mangano 
(« 1 GeV for the 1 5 f g 8 ' > channel), such that the 
averaged cross section is zero below this value. 

Nevertheless, in order to assess the significance 
of the cut-off on the values of the long-distance 
parameters we performed new fits to Tevatron 
data by changing p™ m . We found that, by vary- 
ing p r t nm in the [1,2] GeV range, the 3 s[ S) ME 
showed almost no appreciable change, whereas 
the 1 S^ 8} + 3 P (8) ME changed by a factor two 
at most @ . 

Therefore, we keep in our subsequent analysis 
the value p™ m = 1 GeV as a self-consistent lower 
cut-off. 

3. ALTARELLI-PARISI EVOLUTION 

At high transverse momentum, gluon fragmen- 
tation via the color octet mechanism becomes 
the dominant source of charmonium production. 
On the other hand, AP evolution of the splitting 
gluon produces a depletion of its energy which 
has to be taken into account. If not so, the long- 
distance parameter for the 3 s[ 8 ^ channel would 
be underestimated. 

In this work the AP evolution of the fragment- 
ing gluon was carried out from the evolution of 
the gluonic partner of the J/ ip in the final state 
of the production channel 

9 9 -> 9*{r+JH>) 9 (1) 

the technical reason being that the splitting gluon 
g* actually is not generated in our code ||]. (No- 
tice however that the MEs used in the gener- 
ation include the long-distance evolution of the 
(cc) bound state ||.) 

Indeed, on the average the fragmenting gluon 
should evolve in the transverse plane similarly to 
the other final-state gluon in (1), once corrected 
its virtuality to become of the order of the char- 
monium mass. 

Thereby we get a correcting factor to be applied 
event by event to the transverse momentum of the 
generated J/ip (for the 3 S'j 8 ' ) channel only): 

. AP,* 2 + ml , , 



where Pt (P t * ) is the transverse momentum of the 
final-state gluon without (with) AP evolution. 
At high p t , 

P? P = x p x p t (3) 

where p t is the transverse momentum of the J / ip 
as generated by PYTHIA (i.e. without AP evo- 
lution) . 

Although the above way to implement AP 
evolution may be somewhat indirect, it remains 
in the spirit of our whole analysis, i.e. using 
PYTHIA algorithms whenever possible. In fact 
it provides an energy depletion of the fragment- 
ing gluon in accordance with Cho and Leibovich's 
work |||. 
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Figure 1. Curves obtained from PYTHIA for J/ip 
production at the Tevatron, including AP evolu- 
tion and initial-state radiation on for the CTEQ 
2L parton distribution function. The charm mass 
was set equal to 1.48 GeV. (i) dotted line: Color- 
singlet model, (ii) dashed line: 1 S^ 8) + 3 P (8) , (iii) 
dot-dashed line: 3 5' 1 8 ' ) , (iv) solid line: all contri- 
butions. 



Figure 2. Curves obtained from PYTHIA for 
J/ip production at the LHC, without AP evolu- 
tion and initial-state radiation on for the CTEQ 
2L parton distribution function, (i) dotted line: 
Color-singlet model, (ii) dashed line: 1 S , q 8 ' ) + 3 
P (8) , (iii) dot-dashed line: 3 sf\ (iv) solid line: 
all contributions. 



Figure 3. Curves obtained from PYTHIA for 
J I if) production at the LHC, including AP evolu- 
tion and initial-state radiation on for the CTEQ 
2L parton distribution function, (i) dotted line: 
Color-singlet model, (ii) dashed line: + 3 
F (8) , (iii) dot-dashed line: 3 S[ S) , (iv) solid line: 
all contributions. 



4. NEW EXTRACTION OF COLOR- 
OCTET MEs 

We have performed new fits to the Tevatron 
data M using three PDFs (as in ||) but now 
with AP evolution implemented in the generation 
as described in the preceding section. In table 1 
we show the numerical values for < 8 ( 3 S'i) > 
and the linear combination M 3 = 3 x ( <°s( 3po)> + 
<o 8 (^5 )> ^ o^ained in our analysis. The new 
values for the < 0${ 3 Si) > matrix elements have 
been increased by a factor about 3 w.r.t. AP 
off j|. Conversely the M3 value is slightly de- 
creased, moreover remaining quite smaller than 
in other analysis || . Therefore we conclude that 
the M3 value obtained in our work is reliable and 
in a better agreement with HERA results on J/ ip 
photoproduction ]L3| . 



5. J ftp DIRECT PRODUCTION RATE 
AT THE LHC 

In order to get an estimate of the expected 
charmonium production rate at the LHC, we ran 
PYTHIA with the color-octet model implemented 
in. In figure 3 we show the predicted differential 
cross section times the muonic branching fraction 
of J/ip with AP evolution included. In figure 2 we 
show the same curve with AP off. By comparison, 
one may conclude that, once adjusted the long- 
distance parameters to fit Tevatron data, there 
is no significant difference between both genera- 
tions. 

Let us stress that because of theoretical un- 
certainties associated to different choices for the 
charm mass, energy scales, parton distribution 
function the theoretical curves of figures 2 and 
3 have to be considered as order-of-magnitude 



Table 1 

Color-octet matrix elements (in units of 10~ 3 GeV 3 ) from the best fit to Tevatron data on direct J/ip 
production for different PDFs. Initial-state radiation of gluons and AP evolution as explained in the 
text were swichtcd on. For comparison we quote the values given by Cho and Leibovich: (6.6±2.I) and 
(66±15) respectively. 



ME: 


<0 8 ( 3 5!) > 


M 3 = 3x( <0 <?°> + 


<o 8 Cs )>\ 

3 ) 


CTEQ2L 


9.6±1.5 


13.2±2.1 




MRSDO 


6.8±1.6 


13.2±2.I 




GRVHO 


9.2±1.1 


4.5±0.9 





predictions. Nevertheless the relatively high pro- 
duction rate at high-pt (of the order of the pb 
at p t = 50 GeV) makes especially interesting the 
analysis of charmonia production at the LHC. A 
preliminary study on Y(IS') prompt production 
leads to a similar conclusion for the bottomonia 
family. 

6. CONCLUSIONS 

After checking our analysis by varying the 
lower p t cut-off set in the Monte Carlo generation 
to avoid the problematic p t — >0 region, we con- 
clude that our former choice p™ m = 1 GeV was 
justified and our study of charmonium hadropro- 
duction trustworthy. On the other hand, once im- 
plemented Altarelli-Parisi evolution in our frame- 
work, the new extracted values of the color-octet 
ME < Og( 3 Si) > increase (table 1), recovering a 
similar result as in previous extractions On 
the other hand, the numerical values obtained for 
M3 even decrease slightly, reinforcing the conclu- 
sions reported in our former work (lJ-Q. 

With respect to the extrapolation of the J/tp 
production rate up to LHC energy, the theoret- 
ical prediction incorporating AP evolution (fig- 
ure 3) does not differ significantly from the curve 
obtained without AP evolution (figure 2). This 
means that, from a practical point of view, there 
is no need to modify the code for a fast gener- 
ation as shown in || if the corresponding color- 
octet matrix elements are duly employed. Finally 
we conclude that the analysis of charmonia (and 
bottomonia) prompt production at the LHC de- 
serves special attention in its own right. 



Acknowledgments 

I thank the LHC Workshop b-production sub- 
group for interesting discussions and suggestions. 

REFERENCES 

1. M.A. Sanchis-Lozano and B. Cano-Coloma, 
Nucl. Phys. B (Proc. Suppl.) 55A (1997) 277. 

2. B. Cano-Coloma and M.A. Sanchis-Lozano, 
Phys. Lett. B 406 (1997) 232. 

3. B. Cano-Coloma and M.A. Sanchis-Lozano, 
Nucl. Phys. B 508 (1997) 753. 

4. M.A. Sanchis-Lozano, Nucl. Phys. B (Proc. 
Suppl.) 75B (1999) 191. 

5. E. Braaten and S. Fleming, Phys. Rev. Lett. 
74 (1995) 3327. 

6. T. Sjostrand, Comp. Phys. Comm. 82 (1994) 
74. 

7. G.T. Bodwin, E. Braaten, G.P. Lepage, Phys. 
Rev. D 51 (1995) 1125. 

8. CDF Collaboration, F. Abe at al., Phys. Rev. 
Lett. 69 (1992) 3704; 71 (1993) 2537; 75 

(1995) 1451 ; 79 (1997) 578. 

9. P. Cho and A.K. Leibovich, Phys. Rev. D 53 

(1996) 6203. 

10. M. Cacciari and M. Kramer, Phys. Rev. Lett. 
76 (1996) 4128. 

11. M. Mangano, private communication. 

12. M.A. Sanchis-Lozano, 2nd meeting of the b- 
production working group, LHC99 Workshop, 
April 1999. 

13. B.A. Kniehl and G. Kramer, Eur. Phys. J. C 
6 (1999) 493. 



